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Abstract 

The thermal history of a large class of running vacuum models in which the effective cos¬ 
mological term is described by a truncated power series of the Hubble rate, whose dominant 
term is A(i7) oc 77"+^, is discussed in detail. Specifically, by assuming that the ultrarelativistic 
particles produced by the vacuum decay emerge into space-time in such a way that its energy 
density pr oc T^, the temperature evolution law and the increasing entropy function are an¬ 
alytically calculated. For the whole class of vacuum models explored here we find that the 
primeval value of the comoving radiation entropy density (associated to effectively massless 
particles) starts from zero and evolves extremely fast until reaching a maximum near the end 
of the vacuum decay phase, where it saturates. The late time conservation of the radiation en¬ 
tropy during the adiabatic FRW phase also guarantees that the whole class of running vacuum 
models predicts the same correct value of the present day entropy, Sq ^ (in natural 

units), independently of the initial conditions. In addition, by assuming Gibbons-Hawking 
temperature as an initial condition, we find that the ratio between the late time and primor¬ 
dial vacuum energy densities is in agreement with naive estimates from quantum field theory, 
namely, pao/pai ~ 10“^^®. Such results are independent on the power n and suggests that the 
observed Universe may evolve smoothly between two extreme, unstable, nonsingular de Sitter 
phases. 
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1 Introduction 


A non-singular early de Sitter phase driven by a decaying vacuum energy density was phenomeno¬ 
logically proposed long ago to solve some problems of the Big-Bang cosmology mm- The basic 
idea is closely related to early attempts aimed at solving (or at least alleviating) some cosmic 
mysteries, such as the “graceful exit” problem, which plague many inflationary scenarios (for re¬ 
cent reviews see mm). and also the CCP or cosmological constant problem laEiiTiEiEiiinKn], 
probably, the deepest conundrum of all inflationary theories describing the very early Universe. 

Nevertheless, new theoretical developments are suggesting a possible way to circumvent such 
problems. Results based on the renormalization group (RG) theoretical techniques of quantum 
field theory (QFT) in curved spacetimes combined with some phenomenological insights provided a 
set of dynamical A(H)-models (or running vacuum cosmologies) described by an even power series 
of the Hubble rate [laiisiiii] (cf. mm for a review). In this line, we have discussed in a series 
of recent works a unified class of models accounting for a complete cosmological history evolving 
between two extreme (primeval and late time) de Sitter phases whose spacetime dynamics is 
supported by a dynamical decaying (or running) vacuum energy density [151 13 EH EH E3 [23 E] • 
In such models, the effective vacuum energy density is a truncated power series of the Hubble 
rate, whose dominant term is A(ff) oc where for the sake of generality the power n > 0. 

Unlike several inflationary variants endowed with a preadiabatic phase, this decaying vacuum is 
responsible for an increasing entropy evolution since the very early Universe, described by the 
primeval nonsingular de Sitter spacetime. 

Several theoretical and observational properties of this large class of nonsingular running vac- 
cum scenarios have been discussed in the literature. In particular, Mimoso and Pavon shown their 
thermodynamic consistency based on the generalized second law of thermodynamics (GSLT) by 
taking into account quantum corrections to the Bekenstein-Hawking entropy [22j . Many details 
regarding the late-time dynamics can be found in Refs.[23l |25l |26l [27], and especially in the 
recent, updated and very comprehensive works [28] and [29|. More recently, even the solution of 
the coincidence problem has been detailed discussed in the present framework m. as well as in 
generic decaying vacuum cosmologies m- 

Finally, let us mention that the running vacuum models under study have recently been tested 
against the wealth of SNIa-|-BAO-|-H(z)-|-LSS-|-BBN-|-GMB data - see |32j for a short review - and 
they turn out to provide a quality fit that is significantly better than the ACDM. This fact has 
become most evident in the recent works [331IM] , where it is shown that the significance of the fit 
improvement is at ~ 4(j c.l. Therefore, there is plenty of motivation for further investigating these 
running vacuum models from different perspectives, with the hope of finding possible connections 
with fundamental aspects of the cosmic evolution 

In the present work, we focus our attention on the entropy of the cosmic microwave background 
radiation (GMBR) generated by this large class of non-singular decaying vacuum cosmology. By 
considering that the decaying vacuum process occurs under adiabatic conditions [in the sense that 
the specific entropy (per particle) is preserved] this means that the radiation produced satisfies 
the standard scaling laws, namely, oc and pr oc [36] • Under these conditions, the final 
value of the entropy produced by the decaying vacuum supporting the unstable primeval de Sitter 
phase is exactly the present radiation entropy existing within the current Hubble radius. Within 
this framework, the early decaying vacuum process is not plagued with “graceful exit” problem 
of most inflationary variants and generates the correct number, Sq ~ 10®® [35], regardless of the 
power n present in the phenomenological decaying A (H)-term. In addition, the ratio between the 
primeval and the present day vacuum energy densities is pAo/pAi — 10“^^®, as required by some 
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naive estimates from quantum field theory. 


The article is structured in the following manner. In section 2 we justify the phenomenological 
decaying vacuum law adopted in the paper, whereas in section 3 we set up the basic set of equations 
and the transition from the early de Sitter to the radiation phase is addressed. How inflation ends 
and the temperature evolution law are presented in section 4, while in section 5 we discuss the 
entropy production generated by the decaying vacuum medium. Finally, the main conclusions are 
summarized in section 6. 

2 General model for a complete cosmic history 

The general A(//)-scenario accounting for a complete chronology of the universe (from de Sitter 
to de Sitter) is based on the following expression for the dynamical cosmological term defining the 
relevant class of running vacuum models under consideration [ISKIIKIT]; 

Tjn-\-2 

SttGpa = A{H) = Co + + 3 q . (2.1) 

Here H = a/o is the Hubble rate, a = a{t) is the scale factor, and the over-dot denotes derivative 
with respect to the cosmic time t. By definition pa{H) = A(//)/87r G is the corresponding vacuum 
energy density (G being Newton’s constant). The even powers of H (therefore n = 2,4,...) are 
thought to be of more fundamental origin due the general covariance of the effective action, as 
required by the QFT treatment in curved spacetime [121IT^ ITTl lllj . In the numerical analysis, 
however, we will explore also the cases n = 1, 3 and 4 for comparison. 

The dimensionless free parameters a and v have distinct status. The first can be absorbed (for 
each value of n in the arbitrary scale Hj so that it can be fixed to unity without loss of generality 
(if the scale of inflation is not precisely known) |2r)l I30j . whereas v has been determined from 
observations based on a joint likelihood analysis involving SNe la, Baryonic Acoustic Oscillations 
(BAO) and Cosmic Background Radiation (CMB) data, with the result \v\ = 0(10“^) |23l 1241 [25l 
ESI EH ESI EH] - see especially the most recent analyses in which the = 0 result (associated to 
the ACDM in the post inflationary time) is excluded at ~ 4 (t c.l. [331134j . The small value of v 
is natural since at late times, the dynamical model of the vacuum energy cannot depart too much 
from ACDM. In this connection, by using a generic Grand Unified Theory (GUT), it has been 
shown that \v\ ~ I0“® — 10“^ [13|. Finally, the constant cq with the same dimension of A yields 
the dominant term at very low energies, when H approaches the measured value Hq (from now on 
the index “0” denotes the present day values of the quantities). 

3 Basic equations of the A{H) model 

It is well known that the Einstein field equations (EFE) are valid either for a strictly constant A 
or a dynamical one [mini. Therefore, using the vacuum energy density p\ = A/(87rG) and the 
nominal pressure law p\ = —p\ one can write the EFE in the framework of a spatially flat FLRW 
(Friedmann-Lemaitre-Robertson-Walker) space-time 

R^lu - = 87rGff,^ (3.1) 

where R is the Ricci scalar and T^i, = {pm + Pm)u^.Uu — {Pm — PA)gfiu is the total energy-momentum 
tensor and the index m refers to the dominant fluid component (nonrelativistic matter or radiation). 
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Obviously, in our case, the only difference with respect to the more conventional field equations is 
the fact that A = A{H). In this framework, the local energy-conservation law = 0 which 

insures the covariance of the theory reads: 


PK + Pm + 3(1 + Lo)pmH — 0 , (3-2) 

where we have used Pm = i^Pm for the ordinary cosmic fluid, namely a; = 0 for dust and w = 1/3 
for radiation. In this enlarged framework, the Friedmann equations are given by 

SnGpT = SnCpm + A{H) = 3H^ , (3.3) 

SttGpt = SnGPm - A{H) = -2H - 3H‘^ . (3.4) 


By combining the above equations with the class of vacuum models one obtains the equation 
driving the evolution of the Hubble parameter: 


Co 










(3.5) 


A solution of this equation in the high energy regime [where the term cq/H"^ 1 of (j3.5p can be 
neglected] is given by: 

H(a) = 


[l -I- A) 

where Hj = 

The combination of the EFE and the expression for A yields: 


(3.6) 


PA{a) = PI 
Pr(C') — Pi 


l + V D a3n(l-|-a;)(l-i/)/2 

[l + D a?n{l+u)(l-v)/2-]^ 

(1 — z^)/:) 


p^[a) = PI 


[l -I- ZD 

1 


l+2/n 


\1ln 


(3.7) 

(3.8) 

(3.9) 


The quantity pi in the above equations is the critical energy density defining the primeval de Sitter 
stage 

3H] 


PI 


8ttG' 


(3.10) 


3.1 Prom initial de Sitter stage to radiation phase 

Eor uj = 1/3 the energy density for the vacuum and radiation read: 


PA (a) = Pi 


1 -Fi/ZD a2n(i-i^) 

[l + ZDa2Mi-^)]^+^/'' ’ 


(3.1.1) 


Prip) — Pi 


(1 - i/)ZD 

[l + ZDa2n(l-..)]l+2/’^ ’ 


(3.1.2) 


We can see from (|3.1.1I) that the value (I3.10p just provides the vacuum energy density for a —> 0, 
namely /9 a( 0) = p/. As ji^j 1 we also see that pi/pi ~ thereby effectively showing that 

the constant a can be absorbed in the scale ZZ/, as remarked before. Let us also emphasize from 
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the previous formula that for a —>■ 0 we have Pr/PK oc —>■ 0, i.e. the very early universe 

is indeed vacuum dominated with a negligible amount of radiation. In the rest of the paper, we 
neglect the effects proportional to u (which, since it is the coefficient of is not essential for 
the study of the early universe, the epoch where the is fully dominant). Thus, without 

loss of generality, Hj will be hereafter rescaled so that a = 1 and we set = Oin all the formulae. 
Within this framework, we have Hi = Hj and Ea. (l3.10p becomes 


In addition, from Eq. ()3.6p it follows that the scale factor of the universe takes on an exponential 
form a{t) ~ as long as the condition <C 1 is fulfilled. Obviously, this means that 

the universe is initially driven by a pure nonsingular de Sitter vacuum state, and therefore is 
inflating. However, the above mentioned de-Sitter inflationary phase is only ephemeral. Indeed 
it is easy to check that in the post-inflationary regime, i.e. for a ^ Uj (with E)a3n(i+aj) ^ 
we are led to H oc (or a oc for a; = 1/3 [see Eq. (|3.6I) ]. Therefore, the present model 
evolves smoothly from inflation towards the conventional radiation stage, thereby insuring that 
the initial very large amount of vacuum energy density does not preclude the standard picture of 
the primordial big-bang nucleossynthesis. 

On the other hand, using equation (j3.5p one may check that the decelerating parameter, q = 
—d/aH^ = —1 — H/H^, for arbitrary values of the power index n, reads: 


q{H) = 2 




(3.1.4) 


Naturally, the existence of the radiation stage is not enough to identify precisely the end of inflation 
since the deceleration parameter varied from g = — 1 (de Sitter) to g = 1 (radiation) and the 
inflationary period hnished when g = 0, or equivalently, a = 0. We will discuss this point in the 
next section. 


3.2 When exactly Inflation Ends? 


In order to answer this question we hrst combine Eqs. (I3.1.ip and (|3.1.2p so as to obtain the ratio 
of the radiation energy density (cj = 1/3) to the vacuum energy density: 

. (3.2.1) 

PAW 

In principle, inflation must end when both components - the decreasing vacuum energy density 
and the created radiation energy density - contribute alike. Assuming that the scale factor at the 
point of vacuum-radiation “equality” is a = aeq the above expression implies that D = 1. This 
relation enables us to rewrite the Hubble parameter ()3.6p in the following way: 


H{a, n) 


Hi 

1 -|- (a/Ueq)^” 


1/n 


(3.2.2) 


It follows that the value of the Hubble function at the vacuum-radiation equality depends on the 
value of the parameter n and the initial scale Hi: 


Heq — H (fleq) 


Hi 

2l/n- 


(3.2.3) 
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Figure 1: Upper Panel: Evolution of the normalized Hubble parameter E{a,n) = H{a,n)/Hj 
during the inflationary epoch and its transition into the FLRW radiation era. The Hubble pa¬ 
rameter is normalized with respect to Hj, and the scale factor with respect to aeq, the value for 
which PA = pr (see the text). The lines correspond to the following A{H) oc /Hf scenarios 
[see Eq. m, namely n = 1 (dashed), n = 2 (solid), n = 3 (dot-dashed) and n = 4 (dotted). 
Lower Panel: We provide the relative deviation [1 — E[a,n)/E{a,2)\% of the normalized Hubble 
parameter for the n = 1,3,4 vacuum models with respect to n = 2. 


Now, by inserting this value of H^q into Eq. (j3.1.411 we obtain effectively that q = 0 as should 
be expected. Hence, once the arbitrary scale Hj is fixed, the energy scale or the moment for which 
the inflation ends is readily defined. 


As remarked before, the start of the radiation phase in this model is not characterized by the 
canonical radiation value g = 1, as one might have naively expected. Due to the continuous 
energy exchange between vacuum and radiation there is a short period of time in which q goes 
from g = 0 to the standard result g = 1. In the begin of the standard adiabatic radiation regime 
the deceleration parameter is written as 


<l{Hrad) = 2 




(3.2.4) 


In the approach to the the radiation phase one may safely assume that it started when q{Hj.ad) ~ 
0.9999 with the Hubble parameter given by H^ad = Hi/{2 x 10^)^/” [see Eg. (|3. 2. 411 ]. and from 
Eq. (l3.2.2p we obtain Orad/o^eq — (2 x 10^)^/^”. 
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At this point it is appropriate to make the following comments concerning the cosmic history. 
First, in the full radiation era the value of the scale factor of the universe a^ad becomes at least one 
order of magnitude larger than the corresponding value at the vacuum-radiation equality. Actually, 
the total entropy generated by the vacuum decaying process does not depend on the exact value 
of the ratio Urad/oeg- Second, when the radiation epoch is well left behind the Universe goes into 
the cold dark matter dominated era [Einstein-de Sitter, a{t) oc and, after some billion years 
(~ 7Gyrs) it enters the present vacuum dominated phase, in which A ~ A =const - confer [24t I28j 
and [MUM]. 

In the upper panel of Fig.l we provide the evolution of the normalized Hubble parameter 
E{a,n) = H{a,n)/Hi for the following vacuum models n = 1 (dashed line), n = 2 (solid line), 
n = 3 (long dashed line) and n = 4 (dotted line) [see Eq. (l2.ip ]. We observe that for a <C aeq 
the cosmic evolution begins from an unstable inflationary phase [early de Sitter era, H ~ Hi] 
powered by the huge value Hj presumably connected to the scale of a Grand Unihed Theory 
(GUT). Obviously, when the primeval inflationary era is left behind, particularly for a ^ aeq, 
the cosmic evolution enters smoothly in the standard radiation period H oc o“^. Overall, we 
would like to emphasize that the above natural mechanism for graceful exit is universal for the 
whole class of vacuum models which obey the restriction n > 0. Subsequently when the cq/H"^ 
quantity in Eq. (j3.5h starts to dominate over /Hf (where H Hi) the radiation component 
becomes sub-dominant and the matter dominated era appears. This implies that Ea. ()2.ip reduces 
to E{H) = A+3i'[H‘^ —Hq) which generalizes the traditional AGDM model. In this case the vacuum 
at the present time (cosmological constant) becomes A = 3co -|- 3vHq, where Hq is the Hubble 
constant. More details regarding the late-time dynamics can be amply found in Refs. [23 l I24 p I26 p [27] . 
and especially in the recent, updated and very comprehensive works [28] and |29] . 

Finally, by using n = 2 (corresponding to A{H) oc H^/H^ in the early universe) as a fiducial 
model in the large class (|2.ip . we can appreciate in the bottom panel of Fig.l the relative deviation 
of the normalized Hubble parameters E[a, n) with respect to the n = 2 solution E{a, 2). Obviously, 
when we are far away from the epoch of the vacuum-radiation equality the deviations from the 
fiducial E{a, 2) case are extremely small. On the other hand there is a visible deviation from the 
latter around the transition region a aeq- This deviation becomes at the level of ~ —30% and 
~ 4-20% for n = 1 and n > 2 respectively. 


4 Radiation and its Temperature Evolution Law 


Assuming an “adiabatic” decaying vacuum it has been found that the specific entropy of the 
produced massless particles remains constant, despite the fact that the total entropy can be in¬ 
creasing. This implies that the energy and the number density as a function of the temperature 
are given by the the standard expressions, namely: pr oc and oc T^, but the temperature 
does not obey the scaling relation Tr{t) ~ a{t)~^. Such results were firstly derived based on a 
covariant nonequilibrium thermodynamic description IMllM]), and, more recently, using a kinetic 
theoretic approach [38] • In what follows, we discuss the temperature evolution law in the present 
framework along these lines. 

Let us start with Ea. (l3.1.2p which is rewritten in terms of Ueq as follows 


2n 


_ _ (O'lo^eq) _ _ ^ rpA 


(4.1) 


where in the last equality we included the degrees of freedom (d.o.f.) of the created massless 

^More precisely, the constancy of the specific entropy (per particle) of the produced particles defines the “adia¬ 
batic” vacuum decaying process (see Refs. |36l I37p . 
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modes through the 5 *-factor (see [35]). Now, by solving for the temperature we find: 


n+2 

Tr = 2~^Teq 


(a/Oeq)"/^ 

n+2 5 

[1 + [a/aeqY^] 4 - 


(4.2) 


where T^q = Tr{aeq). 
which is dehned by 


Obviously, T^q is the maximum value of the radiation temperature (|4.2I) 


1/4 


T,g=(}^y -f 


\ 2+5*7 VI6+G5* 


1/4 


(4.3) 


As it is expected Teg is given in terms of the arbitrary initial scale of the primeval de Sitter phase 
(pr, or equivalently Hj). Unlike the value of H(,q (see Eg. (13. 2. 31) 1. this value of Egg is valid for the 
whole class of models since is does not depend on the power n. In particular, this unique maximum 
temperature suggests that the total entropy generated within the horizon and presently observed 
is basically the same for all models (see next section). 

In the top panel of Fig. 2 we present the temperature evolution for several values of n. Note 
that in the very beginning of the evolution (when a 1 —>■ 0 ), the temperature of the created photons is 
also zero in accordance with the fact that = 0 (see Eq. gU). However, for finite values of a 7 ^ 0, 
we observe the existence of two regimes. In the first, a <C Ogg, the radiation temperature (14.21) 
increases as oc (it is linear for n = 2 [l9]), reaching of course its maximum value at a = a^q- 
In this non-adiabatic regime, the vacuum instability guides constantly the model to the standard 
radiation epoch from a process that started in the non-singular de Sitter phase. Note that the 
evolution is different from inflationary models where a highly non-adiabatic “reheating” process 
happens immediately after the adiabatic evolution of the inflaton field (see |35j and references 
therein). In the opposite regime, a ^ aeq, we are well within the radiation epoch when all running 
vacuum models decrease in the classical way, that is, following an adiabatic scaling law, oc a~^. 
As we shall discuss bellow, the power-law increasing of the temperature up to the vacuum-radiation 
equality, beyond of which the universe enters in the standard temperature regime, is the reason of 
the large radiation entropy observed in the present epoch. 

For the purpose of the present study it is also important to find a way to calculate the primeval 
parameters from first principles and indeed one may use the following approaches (from 

now on we consider natural units, H = ks = c = 1). A hrst possibility is to consider that the initial 
de Sitter energy density pi is related with the Planck scale Mp through pj = Mp. An alternative 
approach is to connect pi to the GUT energy scale Mx 10^® GeV and thus pi = (see Refs. 
[201 [ 21 ]). A third possibility is to use the event horizon (EH) of the de Sitter space-time. From the 
analysis of Gibbons & Hawking [39] we know that the temperature of the de Sitter EH in natural 
units is Tqh = H//27r, where Hj is the (constant) Hubble parameter at the de Sitter epoch. In 
our case, following the third approach and combining the last equality of Eq. (14.11) and Eq. (I3.1.3h 
we obtain 


Ti = 


El 

27r 


45 

-) Mp , 


135 , ,4 

pj — -—Mp . 

2g* 


(4.4) 


Then inserting the above expressions into Ea. (l4.3p it is easy to show that the characteristic radi¬ 
ation temperature T^q becomes 


T — 

-‘-eq — 




(4.5) 


It should be stressed that all the above characteristic scales are slightly below the corresponding 
Planck scale, which means that we are inside in the semi-classical QFT regime. For instance, by 
taking 5 * = 106.75 which corresponds to the particle content of the standard model of particle 
physics we find from Ea. (l4.5p that Teg ~ 0.26 Mp which as expected does not depend on the power 
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re. In this respect if we take into account the number of light d.o.f in the GUT then we find that 
the characteristic temperature is still smaller (nearly 10% of Planck mass). At this point it should 
be stressed that the primordial Gibbons-Hawking thermal bath was used only as a peculiar initial 
condition to fix the arbitrary scale Hj. As we shall see in the next section, the ratio between the 
very early and late time vacuum energy densities depends only on the pair characterizing 

the extreme de Sitter phases. It has the expected magnitude thereby also contributing to alleviate 
the so-called cosmological constant problem in the context of the such models (see Eq. (j5.7jP . 


5 The Radiation Entropy 

The total entropy of the radiation included in the present Hubble radius {dp — reads: 

So = ^ gs,o T^o ~ 2.3h-3l0®" ~ , (5.1) 

where T^o — 2.725°K ~ 2.35 x 10“^^ GeV is the CMB temperature at the present time and 
Hq = 2 . 133/1 X 10“^^ GeV (with h ~ 0.67) is the present day Hubble parameter 

In order to check the viability of our model we need to take into account that the total entropy 
should be measured from the initial entropy generated by the decaying vacuum. Since in our 
vacuum model the BBN proceeds fully standard, the equilibrium entropy formula remains valid 
because only the temperature law is modihed [Ml Ea EH]. Therefore, the radiation entropy per 
comoving volume is given by the well known expressions: 


Sr 


f Pr Pr\ 3 4 Pr 3 

\ Tr J 3 Tr 


27r^ o o 
-j^PsTra 


(5.2) 


where Ps is the entropy factor at temperature (at very high temperature Ps is essentially equal 
to the effective number of massless species, p*. However, for lower values there is a correction 
related to the freeze out of neutrinos and electron-positron annihilation). 

With the help of equation (14.2p the comoving entropy (15.2p can be expressed as a function of 
the scale factor as follows: 


Sr 


2(7n+6)/An^2g^ 

45 


Teqalqfn{r) , 


(5.3) 


where r = aja^^q and the function /n(r) depends on the parameter re as given below: 


3(ra+2) 

= 3(n+2) • (5-^) 

^ 4 j 2 

The obtained result for the comoving entropy boils down to the one derived in Ref. m for re = 2, 
as it should. Note also that lima_>.o Sr = 0, as it ought to be expected from the fact that the 
initial de Sitter state is supported by a pure vacuum (no radiation fluid). We also see that during 
the inflationary phase (a <C Oeq), that is, at the early stages of the evolution, the total comoving 
radiation entropy of the Universe increases very fast; in fact, proportional to . For instance, 

for re = 2 (corresponding to 77^-driven inflation) the initial entropy raises as ~ a®. Note also that 
for a = tteq, fn{r = l) = and the associated value S{aeq) is still not the total comoving 

entropy that the decaying vacuum is able to generate (see discussion below Eq. (j3.I.4M . This 
occurs only when a = Urad so that = {arad/<7‘eqY^ ^ 1 and fn{r) ^ 1 for all practical purposes. 

^Following standard lines, we have assnmed the coefficient gs,o = 2 -|- 6 x (7/8) (7i/,o/Tro)^ ~ 3.91 is the entropy 
factor for the light d.o.f. at the present epoch in which we have used T/.o/Uo = (4/11)^^® |35| 
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Figure 2: Top Panel: The evolution of the radiation temperature (normalized with respect to 
its maximum value) during the inflationary period for several values of the free parameter n (see 
caption of Fig.l for definitions) and its transition into the FLRW radiation era. As in Fig.l, the 
lines correspond to the following scenarios n = 1 (dashed), n = 2 (solid), n = 3 (dot-dashed) 
and n = 4 (dotted). Bottom Panel: The evolution of the normalized comoving entropy from the 
inflationary period (where it increases) until reaching the saturation plateau for aja^ad > 1- The 
asymptotic value corresponds the total entropy at present. 


At this point the generated comoving entropy Sr reaches the final value, S^^^ = Sr{arad)i when 
the standard radiation phase begins. 

It thus follows that the asymptotic (adiabatic) value of (15.31) . defined by /(r) ~ 1 for r S> 1, is 
given by: 


= 

‘^rad 


2{7n+6)/4n^2g^ 


45 


-^eq^eq 5 


(5.5) 


a saturated value that must be compared with the present day entropy since the subsequent 
evolution of the model is isentropic 

In the bottom panel of Fig. 2 we show the entropy as a function of the ratio r = aja^q for 
several values of n. Notice, that the entropy is scaled to its value at the vacuum-radiation equality. 
Initially, for a <C Ogg the amount of entropy differs among the vacuum models but when a —>• Oeq 
the corresponding entropies start to converge and subsequently they reach a plateau, namely 
Sr(a) / Sr(aeq) —)■ 8(^+2)/4n fpj. q, ^ characterizes the standard adiabatic phase, which is 

sustained until the present days because the bulk of the vacuum energy A(H) oc already 
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decayed. 

Armed with the above expressions we now compute the prediction of the total entropy inside 
the current horizon volume ~ Using the temperature evolution law [see Eq. (j4.2p ] one may 

see that the expression T^qa^q which appears in the final entropy as given by (15.51) is equal to 
2-3(n+2)/4n2i3^^^3^^^ where Tj-ad = Tr{arad) a^d Qrad was defined in Sect. 13.II [see discussion below 
Eq. (j3.2.4p ]. Note that the n-dependence cancels out and we arrive at to the same final result (the 
power n is important only in the inflationary phase since it determines the time scale in which the 
radiation equilibrium phase is attained): 

SL = (^) T?ad 4ad = So , (5.6) 

where Sq is given by (|5.ip . In the last step we used the entropy conservation law of the standard 
adiabatic radiation phase, which implies that gsT^ad4ad ~ 9s,oT^o4- 


It should be stressed that in the very early de Sitter era, the radiation entropy is zero. However, 
it increases steadily as Sr ~ and, finally, as shown in the bottom panel Fig. 2, deep inside 

the radiation stage becomes constant and approaching its asymptotic present day observed value, 

5o. 

In other words the running vacuum model provides an overall past evolution to the present 
ACDM cosmology by connecting smoothly between two extreme cosmic eras (early inflation and 
dark energy) driven by the vacuum medium [see Eq. (l3.5l) ]. Specifically, using Eq. (14.4h we find that 
the ratio between the associated vacuum energy densities becomes 


PAO _ PAp _ Hp _ HqHaq 
PAi PAi H] 


9 * 

ISOvr M|, 


PAO ^ 10 ^‘^^PAI , 


(5.7) 


in agreement with traditional estimates based on quantum field theory (see Refs. [5] and 0)- Note 
that the late time de Sitter scale, Hp = HqQaq, was used in the above expression [30]. To conclude, 
the decaying vacuum model explains the amount of the radiation entropy and simultaneously it 
also alleviates the so-called cosmological constant problem. 

Needless to say, for a final resolution of the problem one needs to understand the ultimate 
origin of the current value of the cosmological constant. Remarkably, the obtained description of 
the cosmic history is based on a unified dynamical A(77)-term accounting for both the vacuum 
energy of the early and of the current universe. An alternative approach to such description, based 
on the Grand Unified Theory framework, can be shown to provide similar results, see [2ni Ei]. 
Somehow this shows that the obtained results are truly robust and independent of the initial 
conditions reigning in the primeval Universe. 


6 Discussion and Conclusions 

In this paper we have addressed a fundamental issue concerning the thermodynamics of the early 
Universe. It is well-known that within the context of the concordance cosmological model, or 
ACDM model, the thermal history is incomplete and it leads to inconsistencies with the present 
observations. One of the main problems is the well-known horizon problem, which can be rephrased 
thermodynamically as the entropy problem. The concordance model, in fact, cannot offer an 
explanation for the large entropy enclosed in our Hubble sphere, which is tantamount to say that 
it cannot explain the large amount of causally disconnected regions contained in it. 

The well-known solution to these problems is inflation, a patch that has to be added to the 
incomplete ACDM description. While in the more traditional approach inflation is accomplished by 
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postulating the existence of a new fundamental scalar field called inflaton, in the present work we 
have proposed an entirely different (but no less efficient) framework. It is based on the properties 
of a large class of non-singular decaying vacuum models whose structure is that of a truncated 
power series of the Hubble rate, A{H). The involved powers to describe inflation must be higher 
than since the latter can be relevant only for the current Universe. 

In this work we have explored an inflationary dynamics where the decaying vacuum is triggered 
by an arbitrary higher power of the Hubble rate, [ji > q) recently proposed [El flHl flT]. 

In such unified model of the vacuum energy, the effective cosmological term is a dynamical 
quantity, which evolves extremely fast in the early Universe and goes through an approximate de 
Sitter phase in our time - the dark energy epoch. In principle, the late time A(fl)-Universe as 
described by Eq. (2.1) remains very close to the concordance model, but it stills carries a mild 
vacuum evolution (hence a mildly evolving cosmological term) compatible with observations. In 
principle, such a term may act as a smoking gun of its lengthy and energetic history; indeed, 
a much richer history than that associated to the idle A-term inherent to the ACDM model. 
However, since the late time entropy production is very small, for the sake of simplicity, we have 
taken the v parameter equal to zero. Therefore, the model discussed here can be seen as a primeval 
nonsingular phase of the standard ACDM model. Once we leave well back the early times, the v- 
parameter recovers an important role which, in point of fact, makes the running vacuum models 
not only compatible with the current cosmological data but fully competitive with the ACDM 
description |33l El] . 

We have studied in detail some important thermodynamical aspects of that class of dynamical 
vacuum models. Most noticeably we have focused on the issue of the radiation entropy, its origin 
and generation in the first stages of the primeval Universe, and then its final impact on the 
current epoch. Our calculations were based on the assumption that the produced radiation from 
vacuum decay satisfy the standard relations, oc and pr oc [36], a hypothesis related with 
the idea of an “adiabatic” decaying vacuum and the fact that the specific entropy is preserved 
during the process [36|. The basic result is that at early times the temperature of the radiation 
increases (T^ oc a”/^) until its maximum value determined by the equilibrium temperature T^q 
of the vacuum-radiation transition (see Fig. 2) and the same happens with radiation energy 
density. As a consequence, the entropy also increases at very early stages {Sr oc being 

later on conserved during the radiation epoch (neglecting the photon entropy produced in the 
electron-positron annihilation). In this context we have found that the large amount of radiation 
entropy now {Sq ~ ]^g87-88 natural units) can be fully accountable in our dynamical vacuum 
context. We have first shown that the entropy was produced during the inflationary process 
itself at the expense of the continuous vacuum decay. Subsequently, its production stagnated and 
this occurred shortly after the vacuum had lost its energetic power and the Universe entered the 
standard radiation phase. From this point onwards the adiabatic evolution of the cosmos carried 
the comoving entropy unscathed until our days. Overall, the wide class of running A(H)-vacuum 
models provides not only an alternative scenario for inflation (beyond the traditional inflationary 
scalar field models, some of them in serious trouble after the analysis of Planck 2015 data [40] ). 
but also a new clue for graceful exit, which is indeed fully guaranteed within the A(H)-cosmology 
context and does not depend on the power n of the Hubble rate. The remarkable independence 
of both the graceful exit and the entropy prediction from the power n singles out such class of 
dynamical vacuum models from the rest. 

Finally, as originally pointed out, the thermodynamical solution insures that no horizon prob¬ 
lem exists because all points of the current Hubble sphere remained causally connected as of 
the early times when the huge entropy was generated by the decaying dynamics of the primeval 
vacuum. Interestingly enough, the above features are not only universal for the entire class of 
A(H)-models but also independent of the initial conditions of the early Universe. They provide 
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a rather robust basis for the dynamical A(i7)-cosmology, which is currently being tested and will 
continue being tested against the increasingly accurate observations. 
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A FRW radiation phase: Wien’s law and the relation Trad/Tj 

In this appendix a simple argument based on the equilibrium Wien’s law is adopted to show 
that the decaying vacuum drives the model progressively to the radiation phase without reheating 
period (no exit problem). 

To begin with we recall that the initial conditions in our picture are = 0 and Tgh = 
Ti = Hj/27r. However, due to the evolution of the Universe, the Hubble parameter and the 
horizon temperature diminish while the temperature of the created radiation increases due to the 
continuous vacuum decay. 

The model evolves out of equilibrium because the entropy generation is concomitant with the 
inflationary process. In principle, the standard radiation FRW phase is reached when the Wien law 
becomes strictly valid. In what follows we show that it provides an useful constraint on the value 
of the temperature in the begin of the FRW phase thereby suggesting the physical consistency of 
the model. 

In natural units, the standard Wien’s law {XmT = 0.290 cm. K) reads: 

XmT = 1.27. (A.l) 

In order to equalize the horizon temperature, the wavelength at the maximum black body 
intensity (in the begin of the FRW phase, H = Hj-ad, T = T^ad) is expected to be Xm < 

Hence, the Wien’s law takes the form; 


1.27 = XmT < Trad/Hrad- (A.2) 

Now, in our model we know that H^-ad H//(2.10^)i/” (see section 3.2) and Hi = 2'kTi. Thus, it 
follows that 

1.27 < {A,3) 

IttIi 

which can be rewritten as: 


Trad 7.98 
Ti ( 2 . 104 ) 1 /^ 


(A.4) 


The above relation shows two interesting aspects of the model: (i) Trad is much smaller than 
Tgq = Til2^l'^., but it can assume values not dramatically too low in comparison with the both 
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characteristics scales (T/ and T^g); a result in agreement with the demonstrated progressive ap¬ 
proach to the FRW phase (see Eq. (|4.2p and the associated comments), and (ii) it also suggests 
that there is no exit problem in our model (the reheating process is not needed for this class of 
decaying vacuum models). For instance, by taking n= 1, 2 we see that Trad > 4.10 “^Tf, and 
Trad > 0.056Tf, respectively. Note also that the inequality is still more safely satished for values 
of n < 2. This is comprehensible because inflation ends faster for higher values of n. 
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